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Abstract 

The orphan receptor GPR30, previously classified as a G protein-coupled estrogen 
receptor (GPER), has been a subject of debate regarding its ligand specificity. Through 
an integrative approach combining structure elucidation, biochemical binding, and cell 
signaling assays, we demonstrate that estrogen does not directly bind to or activate 
GPR30. Cryo-EM structures of GPR30 reveal an unexpected hydrophilic ligand- 
binding pocket, with striking differences from classical hydrophobic steroid-binding 
sites, inconsistent with estrogen binding. We further confirmed hydrophilic agonists 
like Lys05 as true activators of GPR30, providing structural insights into their binding 
mechanism and receptor activation. Our findings necessitate a paradigm shift in 
defining GPR30's role in estrogen signaling, indicating that its activation occurs 
through mechanisms independent of estrogen binding. This study opens new avenues 
for developing targeted GPR30 ligands and reinterpreting its role in estrogen-mediated 


processes. 
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Introduction 

Estrogen is a pivotal hormone in human physiology, orchestrating a wide range of 
biological processes ranging from reproductive functions to cardiovascular health and 
bone integrity!. Its effects are traditionally understood to be mediated through nuclear 
estrogen receptors, ERa* and ERB>®, which function as transcription factors regulating 
gene expression'*. However, the rapid, non-genomic actions of estrogen imply the 
existence of membrane-associated estrogen receptors enabling swift cellular 
responses’”. G protein-coupled estrogen receptor (GPER)!®, also known as GPR30!"" 
\4 was identified as one such candidate over 18 years ago!*'®. Initial studies reported 
that GPR30 binds estrogen'*'® and mediates rapid estrogen signaling through various 


15-18 including downstream pathways mediated by the G protein subtype, Gq!”, 


pathways 
which induces the generation of the second message, inositol phosphate IP1, leading to 
activation of PKC and calcium signaling. This sparked great interest and debate over 


the role of GPR30 in estrogen physiology”””!. 


The discovery of GPR30 as a novel estrogen receptor (GPER) has provided a 
mechanism of the rapid action of estrogen via GPR30 and facilitated the development 
of GPR30 modulators”. While the endogenous ligand 17-estradiol (E2) activates 


22.23 and raloxifene”*, 


GPR30!5!6, selective estrogen receptor antagonists like tamoxifen 
ER-degrading agent fulvestrant!®, have also been reported to display GPR30 agonism. 
Additionally, several synthetic ligands, including G-1°, were reported to be selective 


GPR30 agonists. 


However, subsequent studies have produced inconsistent results about whether GPR30 
truly serves as an estrogen receptor with independent signaling capacity from nuclear 
estrogen receptors!?°?7, A recent unbiased screen found no activation of GPR30 by 
estrogen, 4-hydroxytamoxifen, or the selective GPR30 agonist G-1 in cells 
overexpressing GPR30, confounding GPR30 as a bona fide estrogen receptor. Instead, 


it identified chlorhexidine, Lys05, and 9-aminoacridine as potent GPR30 activators”. 
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Furthermore, knocking down GPR30 in MCF-7 cells expressing estrogen receptors 
(ERs) and GPR30 showed no impact on E2-induced the generation of cyclic AMP 
(cAMP), the second message downstream of Gs activation, but E2-induced cAMP 
production is affected by double knockout of ERa and ERB”. The presence of 
membrane ERs (mERs)’®, albeit at a limited level of about 3-10% of classical nuclear 
ERs, which are reported to co-operate with GPR307°, further complicated the rapid non- 
genomic signaling events mediated by E23?. There is also evidence that G-1 may signal 
through estrogen receptor variants ERa-36 and ERa-46 instead of GPR303!3%41, In 
addition, rapid signaling by estrogen at the plasma membrane is absence in ER negative 
but GPR30 positive cells *°. Clearly, definitive evidence for direct estrogen binding and 


signaling by GPR30 is still lacking and controversial***4. 


In this paper, we combined biochemical, structural, and functional approaches to 
directly test whether GPR30 interacts with and signals in response to estrogen and 
estrogen-related compounds. We determined cryo-EM structures of GPR30 bound to 
Gq in the presence of high concentrations of estrogen and its related ligands, G-1 and 
fulvestrant, however, none of these compounds was found in these GPR30 structures. 
Functional and radioligand binding assays were consistent with the structural data that 
show no direct GPR30-estrogen interaction. In contrast, the hydrophilic ligand, Lys05, 
was found to bind to the orthosteric pocket of GPR30, revealing a unique mechanism 
of GPR30 activation independent of estrogen binding. By definitively demonstrating 
GPR30 is not a membrane estrogen receptor, our integrative study fundamentally shifts 


thinking and necessitates reinterpreting GPR30’s contributions to estrogen biology. 


Cryo-EM structures contradict GPR30 as a direct estrogen receptor 

GPR30 was hypothetically conceived as a G protein-coupled estrogen receptor that 
directly binds estrogen, offering mechanisms explaining non-genomic effects and 
tissue-selectivity of estrogen signaling. This provocative idea has profoundly impacted 


on our comprehension of rapid estrogen biology. Thus, we initially aimed to use 
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structural biology, specifically cryo-EM, a technique successfully applied in numerous 


GPCR-ligand studies“, to visualize estrogen binding directly to GPR30. 


For our cryo-EM studies, we engineered the human GPR30 by appending a 
haemagglutinin (HA) signal peptide, FLAG, and His tags for expression and 
purification, and introduced a BRIL fusion tag at the N-terminus to enhance complex 
stability (Fig. S1). We utilized an engineered Gaq (see Methods), shown effective in 
mimicking Gaq-coupled GPCR activation in various receptors, like growth hormone 
secretagogue receptor (GHSR)**. Co-expression of GPR30 with this engineered Gag, 
along with rat GB1, bovine Gy2, and scFv16 in Hi5 insect cells, allowed us to form 


stable GPR30-Gq complexes in the presence of E2, G1, or fulvestrant (Fig. S1). 


We determined structures of GPR30-Gq complexes with saturated concentrations of E2, 
the selective agonist G-1, and the ER-degrader fulvestrant at 3.1-3.2 A resolution (Fig. 
S1, Fig. S2, Table S1). In all structures, GPR30 adopts an active-like state typical of 
GPCR-G protein coupling. Surprisingly, despite relatively high resolutions of the 
structures, no ligand density was observed in the GPR30 binding pocket (Fig. 1A-B). 
We additionally solved an apo-GPR30-Ggq structure at 2.9 A resolution (Fig. 1C, Fig. 
S2, Table S1). Structural alignments reveal a striking similarity between GPR30 
conformations in the presence or absence of added ligands, with root mean square 
deviation (RMSD) values of just 0.373 A, 0.378 A and 0.368 A when comparing 
structures solved with E2, G-1, or fulvestrant, to the apo-GPR30-Gq complex, 
underscoring the highly analogous states regardless the presence of these compounds 
(Fig. S3A). A comparison of the EM structures with the AlphaFold® structure showed 
overall concordance yet revealed notable local discrepancies in positioning of 
transmembrane helices, particularly having 3-4 A differences in the extracellular 
regions of TM1, TM2, and TM3 (Fig. 1D), indicating that even for recurrence fold such 
as GPCRs, AlphaFold prediction can still deviated greatly from experimental 


structures””, 
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Notably, examining GPR30’s pocket architecture reveals a striking hydrophilic nature 
enriched with negatively charged residues (Fig. 1E-F). This clashes with 
accommodating hydrophobic estrogen ligands (Fig. S3B), which is in contrast with 
related steroid hormone receptors with hydrophobic pockets as indicated in ERa, ERB, 
and several steroid binding GPCRs (Fig. S3C). In addition, the GPR30 pocket is much 
larger than the corresponding pocket for steroid-related ligands (Fig. S3D). This 
spacious GPR30 cavity also seems unsuitable for smaller estrogen binding. Together, 
these inconsistencies argue that estrogen is unlikely the cognate ligand as originally 
hypothesized based on the cryo-EM structural evidence we established. In the next 
sections, we provide further biochemical and functional evidence that estrogen is not a 


direct ligand of GPR30. 


Functional assays challenge GPR30 as a direct estrogen receptor 
Our structural findings challenge the proposed concept of GPR30 as a membrane 
estrogen receptor. We further investigated this through binding experiments and 


signaling analyses. 


First, we performed the competitive radioligand binding assays using [7H]17f-estradiol 
(PH]-E2) (Fig. 2A) with CV1 cells expressing GPR30 or ERa. As anticipated and in 
alignment with existing literature*!, ERa-transfected cells exhibited high-level binding 
to [SH]-E2 that can be specifically competed by addition of unlabeled E2 (Fig. 2B). 
Conversely, [*H]-E2 showed low basal-level binding to GPR30-transfected cells 
(similar to empty cells) in the presence or absence of increasing concentrations of 
unlabeled E2 (Fig. 2B). E2-related ligands such as tamoxifen and raloxifene (Fig. 2A) 
can also effectively compete with [7H]-E2 in binding to ERa but did not affect the 
nonspecific binding of [*H]-E2 on GPR30 (Fig. 2C-D). We recapitulated these results 


with [>H]-E2 binding to the membrane from Hi5 cells overexpressing GPR30 and ERa, 
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and the consistent results further reinforce the lack of direct GPR30-estrogen interaction 


(Fig. S4A). 


In cell-based functional assays, E2 elicited no changes of cAMP or IP1 in GPR30- 
expressing cells at concentrations up to 100 uM, indicating no stimulation of Gs or Gq 
pathways (Fig. 2E, F). E2 did not reduce the cAMP levels promoted by forskolin either, 
indicating no activation of Gi pathway (Fig. 2G) In contrast, the recently identified 
GPR30 agonist Lys05 (Fig. 2A) induced dose-dependent accumulation of IP1 without 
affecting cAMP levels (Fig. 2E-G), consistent with the knowledge that GPR30 
dominantly couples to Gq!®. Lys05 also induced calcium flux in GPR30 expressing 
cells with an ECso value of ~462 nM, a response not observed with E2-related 
compounds or in parental HEK293 cells (Fig. S4B-C), consistent with activation of the 
Gq-IP1-calctum pathway. Together, our integrated binding and functional data clear 


indicate that GPR30 is not an estrogen receptor. 


Structure of GPR30 bound to Lys05 

Our functional assays have established that Lys05, a molecule characterized by multiple 
positively charged nitrogen atoms (Fig. 2A), effectively activates GPR30, initiating 
downstream signaling through the Gq pathway (Fig. 2E, Fig. S4B-C). Intrigued by this 
interaction, particularly given GPR30's negatively charged ligand-binding pocket (Fig. 


1E), we embarked on a detailed structural investigation. 


We successfully determined the high-resolution structure of the Lys05-GPR30-Gq 
complex at an overall resolution of 2.6 A (Fig. 3A, Fig. S5). This resolution enabled us 
to resolve the entire structure, including the transmembrane helices, intracellular loops, 
and the Gq heterotrimer. Most notably, we observed a distinct EM density for Lys05 
within the GPR30 binding site (Fig. 3A), in stark contrast to our previous GPR30 


structures solved in the presence of E2 and related compounds (Fig. 1 A-B). 
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In the Lys05-bound state, the ligand is nestled within a pocket formed by TM2, TM3, 
TM6, and TM7. This pocket is enriched with negatively charged residues, aligning 
perfectly with Lys05's structure (Fig. 3B), which features two hydrophobic ends linked 
by a positively charged segment (Fig. 2A). This spatial arrangement highlights the 
specificity of Lys05 binding to GPR30. 


A comparative analysis between the apo-GPR30 and Lys05-activated GPR30 structures 
revealed a similar activated conformation (Fig. SSG). The RMSD values for the entire 
complex were 1.223 A and 0.326 A for the receptor alone, indicating a high degree of 
structural similarity. Upon closer inspection, we noted that TM3 and TMS in the Lys05- 
bound structure underwent inward shifts compared to their positions in the apo structure. 
Specifically, TM3 shifted by approximately 0.9 A (Ca of T131°**°) and TMS by around 


2.5 A (Ca of E213°**’), suggesting ligand-induced conformational changes (Fig. S5G). 


Structural analysis reveals a unique distribution of the transmembrane region of GPR30. 
Notably, TM1 and TM4 arrangements in GPR30’s 7TM bundle deviate from 
amine/steroid receptors, such as serotonin receptor, 5HT14° and G protein-coupled bile 
acid receptor, GPBAR™ (Fig. S6A), but match better with peptide-activated GPCRs 
like GHSR* and neuromedin U receptor 2 (NMUR2)™ (Fig. S6B). This architecture 
likely accommodates the Lys05’s linear amine linker. Overall, the Lys05-bound GPR30 
structure provides molecular insights into ligand interactions enabling pharmacological 


targeting of GPR30’s distinctive ligand binding pocket. 


Molecular Recognition of Lys05 by GPR30 

The relatively high-resolution structure of Lys05 coupled with an activated GPR30 
conformation has elucidated the ligand recognition mechanisms in detail (Fig. 4A). 
For analysis, we designated the nitrogen atoms in Lys05's symmetrical structure as N1 


through NS, from the bottom to the top of the ligand structure. Lys05 inserts into the 
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orthosteric pocket lined by TM2, TM3, TM6, and TM7, engaging with the receptor 


through a tri-layered interface (Fig. 4A). 


At the bottom of the binding pocket, the insertion of Lys05 rearranges L137°*? and 
M141°°°. Specifically, L137°°? flips towards TM2, while M141°*°° rotates 
approximately 50° away (Fig. S7A). This rearrangement gives rise to a hydrophobic 
sub-binding pocket consisting of L108**°?, L137°°7, M1413936, W272%48, A3 137”, and 
F314” (Fig. 4A-B). To underscore the importance of these hydrophobic contacts, we 
introduced alanine mutations at these sites and assessed their impact on Lys05 activity. 
Strikingly, alanine substitution at any of these residues, including L1087*7A, L137°*”, 
M1413336, W272°8, and F3147*43, nearly abolished activation by Lys05 (Fig. 4C, Fig. 
S7), without affecting GPR30 surface expression (Table. S2). These data suggest the 


crucial role of the hydrophobic sub-pocket in stabilizing Lys05. 


At the center of the pocket is a hydrophilic layer, where the positively charged N2 and 
N3 amines of Lys05 form polar contacts with E115**°, Q138°°3, E275! and N310”°? 
(Fig. 4A, 4D). Specifically, E1157*© forms polar interactions with the N2 amine of 
Lys05, while N3107*? and E275°*! form polar interactions with the N3 amine of Lys05 
(Fig. 4D). Critically, the N2 and N3 amines provide essential positive charges for the 
accommodation of the surrounding negatively charged pocket (Fig. 4A). As expected, 
E1157°A and Q138°*%A mutations eliminated activation, while E275°°!A reduced 


potency by approximately ten-fold (Fig. 4C, Fig. S7, Table. S2). 


The top of the pocket is a hydrophobic layer that involves weaker interactions between 
the top tryptamine moiety of Lys05 and residues from TM6 and ECL3 of GPR30, 
engaging H300°*'3, H282%58, and P3037? (Fig. 4E). Mutations in these amino acids 
mildly impacted activation (Fig. 4C, Fig. S7, Table. S2), suggesting their auxiliary role 


in LysO5 recognition. 


244 


245 


246 


247 


248 


249 


250 


251 


252 


253 


254 


255 


256 


257 


258 


259 


260 


261 


262 


263 


264 


265 


266 


267 


268 


269 


270 


271 


272 


Activation mechanisms of GPR30 

In the cryo-EM map of the apo-GPR30 Gq structure, several structured waters clearly 
form hydrogen-bonding interactions within the orthosteric pocket. Specifically, two 
well organized water molecules, denoted as W1 and W2, are solved to a local resolution 
of approximately 3.0 A (Fig. S2K), are hypothesized to potentially play a role in 
GPR30's self-activation. These water molecules, positioned in proximity to the 
negatively charged pocket core, actively engage in hydrogen bonds with key residues 
(Fig. 5A-B). W1's positioning and interactions, particularly with E115?*°°, E275°°, 
and N310”°?, coincide with the tryptamine group of Lys05 (Fig. 5C). Meanwhile, W2 
establishes hydrogen bonds with E275°°! and E218°* (Fig. SA-C). These water- 
mediated contacts appear to mimic ligand interactions, promoting GPR30's adoption of 
an active-like state, offering insights into potential basal activation mechanisms of 


GPR30. 


In the absence of an inactive GPR30 structure, we conducted a comparison with the 
inactive K-opioid receptor (KOR)°°, which displays one of the highest sequence 
similarities among all GPCRs, to GPR30. This comparison between active and inactive 
GPCR states illuminated the activation mechanism of GPR30. Notably, the active 
GPR30 structures demonstrate typical rearrangements characteristic of activated 
GPCRs, such as the agonist-induced outward movement of TM6 at the cytoplasmic side, 


similar to changes observed in KOR” (Fig. 5D). 


In the Lys05-GPR30-Gq complex, the binding of Lys05 exerts pressure on the W272*® 
toggle switch, triggering a cascade of conformational shifts (Fig. 5E). This interaction 
induces alterations in the F268°4-V145>° of PIF motif and R155°°° of the DRY 
motif, ultimately disrupting the ionic lock between TM3 and TM6. As a result, TM6 
moves outward by approximately 6.2 A at the R254°°° Ca position. Concurrently, 
TM7's cytoplasmic portion shifts inward by about 2.7 A (Fig. 5D). This inward 


movement allows Y324’°? of NPxxY motif residues, and Y234°"*, to form hydrogen 
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bonds with R155°*°° of DRY (Fig. SF-H). These coordinated structural changes 
facilitate the opening of the cytoplasmic pocket, priming GPR30 for effective coupling 


of downstream G protein. 


Gq coupling of GPR30 

Structural comparisons of the GPR30-Gq complex with other Gq/11-coupled class A 
GPCRs*°** reveal a broadly parallel coupling mechanism. However, our analysis 
identified a distinctive feature in GPR30's coupling with Gaq, particularly in the 
positioning of the aH5 helix of the Gaq protein (Fig. 6A). In the GPR30-Gq complex, 
the N-terminus of «HS is positioned closer to the junction of TM5 and TM6 in GPR30, 
diverging by approximately 15° compared to other Gq-coupled receptors such as 
GSHR*, NMUR2™ and CCKAR* (Fig. 6A). This unique orientation may reflect the 
specific activation characteristics of GPR30, which induces a comparatively smaller 


displacement of TM6 (Fig. 6A). 


The interaction between GPR30 and Gaq involves two critical regions of the receptor: 
the engagement of the Gaq a5 helix by TM3, TMS, TM7, and H8 of GPR30, and 
interactions mediated by GPR30's intracellular loops (ICLs) (Fig. 6B). The major 
interface is formed by the distal aH5 region of Gaq and the intracellular receptor cavity. 
Key hydrophobic and polar interactions at this interface include:1) a hydrophobic stack 
between V394SH5-26, 13936-4525, 1 3886-H5-20 7 384G-H5.16. 13835-45-15 on Gaq a5 helix, 
and GPR30 residues L15994, L241°*°, A356%%, M260, 126137; 2) polar 
interactions such as the Q385°%5!7 - R254°° hydrogen bond and N392°4>.4 
interacting with T330*” in helix 8 and backbone of $3257*™ (Fig. 6B-C). 


Significant ICL2-Gaq interactions are also observed, including the insertion of 
M163*4**! into a hydrophobic cavity on Gag. This interaction is reminiscent of those 
seen in other Gq-GPCRs*$5457 (Fig. 6B and 6D). Additionally, R169457 forms a polar 


QG-H5.22 


interaction with E39 , with its backbone carboxyl group interacting with R38 on 
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aN of Gaq. R164°**? establishes a hydrogen bond with the D215 backbone of Gaq. 


A unique aspect of the GPR30-Ggq structure is the increased interaction between TM6- 
ICL3 of GPR30 and Gag, facilitated by a smaller outward movement of TM6 in GPR30 
(Fig. 6A). Notably, R254°°° on GPR30 forms a stabilizing interaction with Q385°%:!7 
and D381°>!3 on Gaq aH5. Furthermore, R248™'? forms a salt bridge with 
D3789¥5-10 on aH5 (Fig. 6E). These additional contacts, unique to the GPR30-Gq 
complex, highlight the distinct conformational adaptations of GPR30 in its coupling 


with Gaq. 


Discussion 

The study of estrogen signaling is pivotal for understanding a wide range of 
physiological processes, from reproductive health to neurological functions'*. GPR30, 
previously conceptualized as a membrane estrogen receptor (GPER)!”, has been central 


to discussions about rapid, non-genomic estrogen responses?””! 


. However, our 
comprehensive structural, biochemical, and functional analyses necessitate a 
fundamental shift in understanding GPR30's role. We have clearly demonstrated that 
GPR30 does not directly bind or get activated by estrogen, thereby challenging its 


classification as a direct membrane estrogen receptor. 


The cryo-EM structures of GPR30 revealed no density for E2, G-1, or the ER-degrader 
fulvestrant. Instead, the orthosteric pocket of GPR30 exhibits a hydrophilic nature, 
unsuitable for binding of hydrophobic steroid hormones like estrogen. This finding 
aligns with our radioligand binding assays, which also indicates no specific estrogen 
binding to GPR30. Consequently, we must reinterpret GPR30’s role in estrogen 
signaling, considering indirect mechanisms or completely estrogen-independent 


functions of this receptor. 


The definitive exclusion of GPR30 as a membrane estrogen receptor redirects our focus 
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to other potential mediators of rapid estrogen responses. Candidates like ERa-36*°, 
ERa-46*! and membrane-associated ERs (mERs) warrant further investigation to 
understand their roles and interactions in membrane-initiated signaling?’?!*°, 


Clarifying the crosstalk or cooperation between these receptors and GPR30 could 


provide deep insights into estrogen's complex signaling pathways. 


The nature of GPR30's native ligand remains an open question, especially since 
estrogen is now excluded. The hydrophilic pocket architecture of GPR30 and the unique 
arrangement of its 7TM bundle suggest an affinity for ligands differing from traditional 
hydrophobic steroids. We have successfully characterized Lys05, a synthetic small 
molecule, which activates GPR30 and the downstream Gq signaling pathway. The 
extensive interactions between Lys05 and GPR30, as revealed by our structural analysis, 


advance our understanding of GPR30’s activation and signaling mechanisms. 


Our findings have significant implications for the field of hormone signaling and the 
study of GPCRs. The discovery that GPR30 does not act as a direct estrogen receptor 
calls for a reconsideration of its nomenclature. Referring to GPR30 as GPER may no 
longer be appropriate, considering its distinct ligand-binding and activation properties. 
This study lays the groundwork for developing new tools to probe GPR30's function, 
potentially leading to a better understanding of its physiological roles and its 
connections to estrogen biology. Unraveling these mechanisms will offer profound 
insights into hormone signaling complexities and pave the way for innovative 


therapeutic approaches targeting GPR30. 
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Methods 

Cell lines 

Spodoptera frugiperda (Sf9, Expression systems) and Trichoplusia ni. (High Five, 
Thermo Fisher) cells were grown in ESF 921 medium at 27°C and 120 rpm. HEK293T 
cells were grown in a humidified 37°C incubator with 5% CO2 using media 
supplemented with 100 I.U./mL penicillin and 100 mg/mL streptomycin (Invitrogen). 
The human cell lines HEK293T were maintained in DMEM (VWR) containing 10% 
fetal bovine serum (FBS, VWR). 


Plasmids construct 

For structure determination of the E2, G1, fulvestrant and apo GPR30-Gq complex, the 
full-length gene sequence of wild-type human GPR30 was cloned into pFastBac vector 
with an N-terminal HA signaling peptide sequence followed by a FLAG tag, a 10xHis 
tag and a BRIL to facilitate expression and purification’. The engineered Gaq 
(hereafter Gaq) construct was generated based on mini-Gs/q713, which carries two 
dominant-negative mutations (corresponding to G203A and A326S) to decrease the 
affinity of nucleotide-binding and facilitate the stability of GaBy complex*’. The N 
terminal 1—18 amino acids and the a-helical domain of the mini-Gs/q71 were 
substituted by the corresponding sequences of the human Gail contributing to binding 
to the antibody fragments scFv16°. All the three G protein complex components, Gaq, 
rat GB1 and bovine Gy2, were cloned into pFastBac vector separately. For structure 
determination of the Lys05 stimulated GPR30-Gq complex, the engineered Gaq was 
fused to the C-terminus of GPR30 with flexible linker between them. GB1 and Gy2 


were cloned into a pFastBac Dual vector. 


Complex expression and purification 
The complexes are all expressed in Hi5 insect cells (Invitrogen). For expression of E2, 
G1 and fulvestrant stimulated GPR30-Gq complex, and apo-GPR30-Gq complex, cell 


cultures were grown in ESF 921 medium (Expression system) to a density of 3x10 per 
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mL with virus preparations for GPR30, Gaq, GB1, Gy2, and scFv16 at the ratio of 
1:1:1:1:1. For expression of Lys05 activated GPR30-Gq complex, cell cultures were 
grown to a density of 3x10° per mL with two separate virus preparations for GPR30- 
Gaq and GB1y2 at a ratio of 1:1.2. The infected cells were cultured at 27°C for 48 h 


before collection by centrifugation and the cell pellets were stored at -80°C. 


For the purification of the E2, G1 and fulvestrant stimulated GPR30-Gq complex, cell 
pellets from 1 L culture were thawed at room temperature and resuspended in low salt 
buffer containing 20 mM HEPES pH 7.4, 100 mM NaCl, 5 mM CaCh, 5 mM MgCh, 
10% glycerol, protease inhibitor cocktail (Thermo Fisher Scientific). The GPR30-Gq 
complexes were formed on membrane in the presence of 1 uM E2 or 1 uM Gl 
(Targetmol.USA) or 1 uM fulvestrant (Sigma), and treated with apyrase (20 mU mL-1, 
NEB), followed by incubation 1 hour at 4°C. Subsequently, 0.5% (w/v) n-dodecyl-B-d- 
maltoside (DDM, Anatrace) and 0.1% (w/v) cholesteryl hemisuccinate (CHS, Anatrace) 
were added to solubilize complexes for 2-3 h at 4°C. Insoluble materials was removed 
by centrifugation at 30,000 g for 30 min and the supernatant was immobilized by batch 
binding to Talon affinity resin. After batch binding, the TALON IMAC resin with 
immobilized protein complex was manually loaded onto a gravity flow column. The 
TALON IMAC resin was washed with 10 column volumes of 20 mM HEPES, pH 7.4, 
100 mM NaCl, 5 mM CaCh, 5 mM MgCh, 30 mM imidazole, 10% glycerol, 0.1% 
LMNG (w/v), 0.02% CHS (w/v), 1 uM E2 or 1 uM G1 or 1 uM fulvestrant and eluted 
with the same buffer plus 250 mM imidazole, 1 uM E2 or 1 uM G1 or 1 uM fulvestrant. 
The mixture was then purified by SEC using a Superdex 200 10/300 GL column (GE 
healthcare) in 20 mM HEPES, pH 7.4, 100 mM NaCl, 0.00075% (w/v) LMNG, 0.00025% 
(w/v) GDN and 1 uM E2 or 1 uM G1 or 1 uM fulvestrant. For the purification of apo- 
GPR30-Gq complex, all steps are same as aforementioned without ligand added. For 
the purification of the Lys05 stimulated GPR30-Gq complex, 500 uM Lys05 
(Targetmol.USA) was added throughout the purification, and the eluted protein was 


incubated with 20 g/mL of Nb35 at 4 °C for another 1 h. The mixture was then purified 
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by SEC using a Superdex 200 10/300 GL column (GE healthcare) in 20 mM HEPES, 
pH 7.4, 100 mM NaCl, 0.00075% (w/v) LMNG, 0.00025% (w/v) GDN and 500 uM 


Lys05. All sample are finally concentrated to ~8 mg/mL for cryo-EM grids preparation 


Cryo-EM grid preparation and data collection 

For cryo-EM grids preparation of the GPR30-Gq complexes, 3 uL of the protein at ~8 
mg/mL were loaded onto a glow-discharged holey carbon grid (Quantifoil Au 300 mesh 
R1.2/1.3), and subsequently were plunge-frozen in liquid ethane using Vitrobot Mark 
IV (Thermo Fischer Scientific). The data of E2, G1 and fulvestrant stimulated GPR30- 
Gq complex were collected at the electron microscopy facility of Shanghai Institute of 
Materia Medica (SIMM) on a Titan Krios at 300 kV using Gatan K3 Summit detector 
with a pixel size of 1.071 A at the Shanghai Institute of Materia Medica, Chinese 
Academy of Sciences. Images were taken at a dose rate of about 8.0 e/A2/s with a 
defocus ranging from -1.0 to -2.0 um using the Serial-EM software®'. The total 
exposure time was 8 s, and 36 frames were recorded per micrograph. A total of 4,773, 
4,840 and 7,980 movies were collected for E2, G1 and fulvestrant stimulated GPR30- 


Gq complex, respectively. 


For apo GPR30-Gq complex, cryo-EM imaging was collected on a Titan Krios at 300 
kV using Gatan K3 Summit detector with a pixel size of 0.824 A at the Shanghai 
Advanced Center for Electron Microscopy, Shanghai Institute of Materia Medica, 
Chinese Academy of Sciences. Images were taken at a dose rate of about 8.0 e/A2/s 
with a defocus ranging from -1.0 to -2.0 um using the EPU software (FEI Eindhoven, 
Netherlands). The total exposure time was 8 s, and 36 frames were recorded per 


micrograph. A total of 6,691 movies were collected. 


For Lys05 activated GPR30-Gq complex, Cryo-EM imaging was collected on a Titan 
Krios equipped with a Falcon 4 direct electron detection device at 300 kV. Images were 


taken with a pixel size of 0.73 A, a defocus ranging from -1.0 to -2.0 um using the EPU 
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software (FEI Eindhoven, Netherlands). We collected a total of 9,417 movies with total 
dose of 50 e/A2/s over 2.5 s exposure on each EER format movie. Each movie was 


divided into 36 frames during motion correction. 


Image processing and map construction 

For the dataset of the E2 stimulated GPR30-Gq complex, image stacks were subjected 
to beam-induced motion correction using MotionCor2. Contrast transfer function 
(CTF) parameters for each non-dose-weighted micrograph were determined by Gctf®. 
Automated particle selection and data processing were performed using RELION-3.0 
beta2. Automated yielded particles were extracted on a binned dataset with a pixel 
size of 2.142 A and were subjected to reference-free 2D classification, producing 
1,199,129 particles with well-defined averages. Further 3D classification produced two 
good subsets showing clear structural features accounting for 527,631 particles. These 
particles were re-extracted into a pixel of 1.071 and subsequently subjected to CTF 
refinement, Bayesian polishing. Further 3D classifications with no alignment produced 
a high-quality subsets accounting for 177,591 particles. 3D refinement and post process 
generated a map with an indicated global resolution of 3.2 A at a Fourier shell 


correlation (FSC) of 0.143. 


For the dataset of the fulvestrant stimulated GPR30-Gq complex, image stacks were 
subjected to beam-induced motion correction using MotionCor2®. CTF parameters for 
each non-dose-weighted micrograph were determined by Gctf®*. Automated yielded 
particles were extracted on a binned dataset with a pixel size of 2.142 A and were 
subjected to reference-free 2D classification, producing 3,410,733 particles with well- 
defined averages. Further 3D classification produced a good subset showing clear 
structural features accounting for 826,321 particles. These particles were re-extracted 
into a pixel of 1.071 and subsequently subjected to 3D classifications with no alignment 


produced a high-quality subsets accounting for 536,691 particles. 3D refinement and 


23 


690 


691 


692 


693 


694 


695 


696 


697 


698 


699 


700 


701 


702 


703 


704 


705 


706 


707 


708 


709 


710 


711 


712 


713 


714 


715 


716 


717 


718 


post process generated a map with an indicated global resolution of 3.2 A at a FSC of 


0.143. 


For the dataset of the G1 stimulated GPR30-Gq complex, image stacks were subjected 
to beam-induced motion correction using MotionCor2™. CTF parameters for each non- 
dose-weighted micrograph were determined by Gctf**. Automated particle selection 
and data processing were performed using RELION-3.0 beta2™. Automated yielded 
and subsequent reference-free 2D classification produced 1,065,934 well-defined 
particles were imported into cryoSPARC.v4®. Further ab-initio model and 
heterogeneous refinement followed by a resolution ladder hetero-refinement produced 
a high-quality subset accounting for 411,141 particles. The Non-uniform refinement 


generated a map with an indicated global resolution of 3.1 A at a FSC of 0.143. 


For the dataset of the apo-GPR30-Gq complex, the single particle analysis was 
performed with cryoSPARC.v4™. Dose-fractionated image stacks were subjected to 
motion correction by MotionCor2. CTF parameters for micrograph were estimated by 
patch CTF estimation. The auto-picked particles were extracted and subjected to 
reference-free 2D classification. The selected 1,072,840 particles were subjected to ab- 
initio model and heterogeneous refinement, produced a high-quality subset accounting 
for 345,373 particles. Further ab-initio model and heterogeneous refinement produced 
a subset of 300,561 particles for Non-uniform refinement and generated a map with an 


indicated global resolution of 2.9 A at a FSC of 0.143. 


For the dataset of the Lys05 activated GPR30-Gq complex, the single particle analysis 
was performed with cryoSPARC.v4®. Dose-fractionated image stacks were subjected 
to motion correction by MotionCor2™. CTF parameters for micrograph were estimated 
by patch CTF estimation. The auto-picked particles were extracted and subjected to 
reference-free 2D classification. The selected 1,457,109 particles were subjected to ab- 


initio model and heterogeneous refinement, produced a high-quality subset accounting 
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for 274,405 particles. Further ab-initio model and heterogeneous refinement produced 
a subset of 209,371 particles for Non-uniform refinement and generated a map with an 


indicated global resolution of 2.6 A at a FSC of 0.143. 


Model building and refinement 

Predicted model of active-state GPR30 receptor from AlphaFold” were used as initial 
model for rebuilding and refinement against the electron microscopy density map. 
UCSF Chimera-1.14°° was used to dock the model into the electron microscopy density 
map, and followed by iterative manual adjustment and rebuilding in COOT-0.9.6% and 
ISOLDE-1.2°%. Then models were further refined and validated in Phenix-1.20 
programs” (Table S1). Structural figures were generated using UCSF Chimera-1.14%, 
ChimeraX-1.27° and PyYMOL-2.0 (https://pymol.org/2/). 


Cell transfection 

The Wild type GPR30 gene was subcloned into pcDNA3.0 vector with the addition of 
an N-terminal HA tag. All the mutations used for functional studies were generated by 

QuickChange PCR and were verified by DNA sequencing. HEK293 cells were obtained 

from ATCC (Manassas, VA, USA) and cultured in DMEM supplemented with 10% (v/v) 
FBS, 100 mg/L penicillin, and 100 mg/L streptomycin in 5% COz2 at 37°C. For transient 

transfection, approximately 2x10° cells were mixed with plasmids in 200 uL 

transfection buffer, and electroporation was carried out with a Scientz-2C 

electroporation apparatus (Scientz Biotech, Ningbo, China). The experiments were 


carried out 24 h after transfection. 


Surface expression analysis 

24 hours after transfection, cells were washed with PBS, fixed with 4% PFA for 15 min, 
and then blocked with 2% BSA for 1 h. Next, cells were incubated with the polyclonal 
anti-HA (Sigma, H6908) overnight at 4°C and then horseradish peroxidase (HRP)- 
conjugated anti-rabbit antibody (CST, 7074S) for 1 h at room temperature. Then cells 


were washed and incubated with 50 uL tetramethylbenzidine (Sigma, T0440) for 30 
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min before the reaction was stopped with 25 uL TMB Substrate stop solution (Beyotime, 
P0215). Absorbance at 450 nm was quantified using a FlexStation III microplate reader 


(Molecular Devices). 


IP1 accumulation assay 

IP1 accumulation was measured using the IP-ONE Gq HTRF kit (Cisbio, 621PAPEJ). 
Briefly, 24 hours after transfection, cells were harvested and resuspended in DMEM 
containing 0.1% BSA at a density of 4 million per milliliter. 5 ml cells were aliquoted 
into 384-well plate and stimulated with 5ml compounds in DMEM (0.1% BSA) 
supplemented with 50 mM LiCl. After incubation at 37 °C for 30 minutes, 5 ml IP1 d2 
reagent and 5ml IP1 Tb cryptate antibody were added. After incubation at room 
temperature for another 1h, HTRF was read on EnVision multiplate reader 
(PerkinElmer). ECso values for each curve were calculated by Prism 8.0 software 


(GraphPad Software). 


Calcium assay 

Cells were seeded in 96-well plates at a density of 3x104 cells/well and cultured 
overnight. The cells were then incubated with 2 umol/L Fluo-4 AM in HBSS (5.4 
mmol/L KCl, 0.3 mmol/L Na2HPOs, 0.4 mmol/L KH2PO,, 4.2 mmol/L NaHCOs, 1.3 
mmol/L CaCl, 0.5 mmol/L MgCl, 0.6 mmol/L MgSOg, 137 mmol/L NaCl, 5.6 mmol/L 
D-glucose and 250 umol/L sulfinpyrazone, pH 7.4) at 37 °C for 40 min. After a 
thorough washing, 50 uL of HBSS was added. After incubation at room temperature 
for 10 min, 25 uL of agonist was dispensed into the well using a FlexStation III 
microplate reader (Molecular Devices), and the intracellular calcium change was 
recorded at an excitation wavelength of 485 nm and an emission wavelength of 525 nm. 
ECs and Emax values for each curve were calculated by Prism 8.0 software (GraphPad 


Software). 


cAMP accumulation assay 
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cAMP accumulation was measured using the cAMP kit (PerkinElmer, TRF0264). 
Briefly, 24 hours after transfection, cells were harvested and resuspended in DMEM 
containing 0.1% BSA at a density of 2x10% cells/mL. Cells were then plated onto 384- 
well assay plates at 1000 cells/5 uL/well. To test the Gs activity of compounds on 
GPR30, another 5 uL buffer containing compounds with 500 uM IBMX were added 
to the cells. To test the Gi activity of compounds on GPR30, 5 ulL buffer containing 
compounds with 1 uM forskolin and 500 uM IBMX were added to the cells. After 
incubation at 37 °C for 30 minutes, intracellular cAMP level was tested by a LANCE 
Ultra cAMP kit and EnVision multiplate reader according to the manufacturer’s 
instructions. ECso values for each curve were calculated by Prism 8.0 software 


(GraphPad Software). 


Ligand Binding assay 

CV1 cells were cultured in DMEM medium with 10% FBS and seeded at a density of 
30,000 cells/well in Isoplate-96 plates (PerkinElmer). Twenty-four hours after 
transfection with the GPR30 or ERa, CV1 cells were washed twice and incubated with 
blocking buffer (DMEM supplemented with 25 mM HEPES and 0.1% (w/v) BSA, pH 
7.4) for 2 h at 37°C. For homogeneous competition binding, radiolabeled [2H]-E2 
(PerkinElmer, 3 nM) and unlabeled compound at seven decreasing concentrations 
(Estradiol, G-1, Tamoxifen, Raloxifene, 10 uM to 1 pM) were added and competitively 
reacted with the cells in blocking buffer at RT for 3 h. Following incubation, cells were 
washed three times with ice-cold PBS and lysed by 50 uL lysis buffer (PBS 
supplemented with 20 mM Tris-HCl, 1% Triton X-100, pH 7.4). The radioactivity was 
subsequently counted (counts per minute, CPM) in a scintillation counter (MicroBeta 
2 Plate Counter, PerkinElmer) using a scintillation cocktail (OptiPhase SuperMix, 
PerkinElmer). Data were analyzed by nonlinear regression using GraphPad Prism 8.0 


software. 
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Hi5 cell binding assays were carried out in filter plates. His membrane homogenates 
expressing ERa or GPR30 (10 ug protein per well) were incubated in membrane 
binding buffer (20 mM HEPES-NaOH and 10 mM EDTA, pH 7.4) with 50 nM [>H]- 
E2 at room temperature for 3 h. Following incubation, the samples were filtered rapidly 
in vacuum through glass fiber filter plates (PerkinElmer). After soaking and rinsing 4 
times with ice-cold PBS, the filters were dried before the addition of 50 uL of 
scintillation cocktail (PerkinElmer) and counted for radioactivity in a scintillation 


counter (PerkinElmer). Data were analyzed by GraphPad Prism 8.0 software. 
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Figure 1. Cryo-EM structures contradict GPR30 as a direct estrogen receptor. 
(A) Cryo-EM maps of GPR30-Gq structures in the presence E2, G-1, and fulvestrant. 
(B) No ligand density observed in the GPR30 binding pocket. 

(C) Cryo-EM density map and cartoon presentation of the apo-GPR30-Gq complex. 
(D) Comparison of the EM GPR30 structure with the AlphaFold predicted GPR30 
structure. 


(E-F) Negatively charged (E) and hydrophilic (F) nature of GPR30’s binding pocket. 
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Figure 2. Functional assays challenge GPR30 as a direct estrogen receptor. 

(A) Chemical formula of E2, tamoxifen, raloxifene and Lys05. 

(B-D) Radioligand competitive binding results of E2 (B), tamoxifen (C), or raloxifene 
(D) to GPR30, ERa and empty cell. Data are mean + S.E.M. from 3 independent 
experiments (n = 3). 

(E-G) Cell-based functional assays to evaluate stimulation of E2, Lys05, or E2+Lys05 
for Gq (E), Gs (F), or Gi (G) pathways. Data are mean + S.E.M. from 3 independent 


experiments (n = 3). 


Figure 3. Structure of GPR30 bound to Lys05. 

(A) Cryo-EM density map and cartoon presentation of the Lys05-GPR30-Gq complex. 
The complex density map is shown at a contour level of 0.15. 

(B) Positively charged Lys05 accommodated in the negatively charged ligand binding 
pocket of GPR30. 
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Figure 4. Molecular Recognition of Lys05 by GPR30. 


(A) Diagram of Lys05 interacting residues, the background indicated two hydrophobic 
(orange) and one hydrophilic (light blue) moieties of Lys05. The residues and dotted 
line in red color represent polar nature or interactions, and that in black represent 
hydrophobic nature or interactions. 

(B, D, E) Detailed interactions of the tri-layered interface between Lys05 and GPR30, 
the bottom hydrophobic interactions (B), the center hydrophilic interactions (D), and 
the top hydrophobic interactions (E). 

(C) Effects of mutating GPR30-Lys05 interacting residues within the ligand binding 
pocket. The pECso is shown in the Bar graph. Data are mean + S.E.M. from 3 


independent experiments (n = 3). 
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Figure 5. Basal activity and ligand-induced activation of GPR30. 

(A-B) Side view (A) and top view (B) of water molecules and surrounding residues in 
the apo-GPR30-Gq structure pocket. 

(C) Superposition of Lys05 and water molecules in the structures of the GPR30-Gq 
complexes. 

(D) Superposition of activated GPR30 (cornflower blue) with active KOR (gray; PDB 
code: 8F7W) and inactive KOR (orange; PDB code: 4DJH). Notable conformational 
changes occur at intracellular ends of TM6 and TM7 upon receptor activation, side view 
(left) and bottom view (right). 


(E) The “switch” W272™** of GPR30 display relative rotameric change when sensing 


Lys05 and further induced conformational change of F268° 4, 
(F-H) The key D-R155**°°-Y motif (F), P-V145°**? — F268° (G) and N-P’**°-xx- 


Y32475? (H) motifs displayed conformational rearrangement in activated GPR30. 
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Fig. 6 Gaq coupling of GPR30. 


(A) Compared to GSHR (PDB code: 7F9Z), CCKAR (PDB code: 7EZM) and NMUR2 


(PDB code: 7W55), the N-terminus of aH5 is positioned closer to the junction of TM5 


and TM6 in GPR30, diverging by approximately 15°. 
(B) Diagram of the GPR30 - Gaq contacts. 


(C) Key interactions between the TMs that forms the core cavity for Gq coupling and 


aH5 of Gaq. 


(D-E) Key interactions between the ICL2 (D) and ICL3 (E) from GPR30and aH5 of 


Gaq. 
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Fig. S1 Purification and structure determination of GPR30-Gq complexes in the 
presence of E2 and fulvestrant. 
(A) Representative size exclusion chromatography (SEC) profiles and SDS-PAGE 
analysis of GPR30-Gq complex. 
(B-C) Representative cryo-EM image (B) and 2D classification averages (C) of GPR30- 


Gq complex in the presence of E2. 

(D-E) Cryo-EM data processing flowcharts (D) and the Fourier shell correlation (FSC) 
curves (E) of GPR30-Gq complex in the presence of E2. The global resolution of the 
final processed density map estimated at the FSC = 0.143 is 3.2 A. 

(F) The density maps of TM1—TM7 and helix 8(H8) of GPR30 in GPR30-Gq complex 
in the presence of E2. 

(G-H) Representative cryo-EM image (G) and 2D classification averages (H) of 
GPR30-Gq complex in the presence of fulvestrant. 

(l-J) Cryo-EM data processing flowcharts (I) and the Fourier shell correlation (FSC) 
curves (J) of GPR30-Gq complex in the presence of fulvestrant. The global resolution 
of the final processed density map estimated at the FSC = 0.143 is 3.2 A. 

(K) The density maps of TM1—TM7 and helix 8(H8) of GPR30 in GPR30-Gq complex 


in the presence of fulvestrant. 
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Fig. S2 Purification and structure determination of GPR30-Gq complexes in the 


presence of G1 and apo- GPR30-Gq complexes. 

(A-B) Representative cryo-EM image (A) and 2D classification averages (B) of 
GPR30-Gq complex in the presence of G1. 

(C-D) Cryo-EM data processing flowcharts (C) and the Fourier shell correlation (FSC) 
curves (D) of GPR30-Gq complex in the presence of G1. The global resolution of the 
final processed density map estimated at the FSC = 0.143 is 3.1 A. 

(E) The density maps of TM1—TM7 and helix 8(H8) of GPR30 in GPR30-Gq complex 
in the presence of G1. 

(F-G) Representative cryo-EM image (F) and 2D classification averages (G) of apo- 
GPR30-Gq complex. 

(H-I) Cryo-EM data processing flowcharts (H) and the Fourier shell correlation (FSC) 
curves (I) of apo-GPR30-Gq complex. The global resolution of the final processed 
density map estimated at the FSC = 0.143 is 2.9 A. 

(J) The density maps of TM1—TM7 and helix 8(H8) of GPR30 in apo-GPR30-Gq 
complex. 

(K) Local resolution for the density of water molecules (W1—W2) in the ligand-binding 
pocket of the apo-GPR30-Ggq structure. 
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Fig. S3 Comparisons of GPR30 structures, and the difference between GPR30’s 
binding pocket and that of related steroid hormone receptors 

(A) Comparisons of GPR30 structures. 

(B) Chemical formula of hydrophobic ligand E2, G1 and fulvestrant. 

(C-D) The hydrophilic nature and larger volume of the GPR30 pocket is in contrast to 


that of the related steroid hormone receptors. 
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Fig. S4 Binding assay results of E2 to GPR30 and validation of Lys05 induced 
GPR30 activation and Gq signaling. 

(A) Binding assay results of [7H]-E2 to GPR30 in Hi5 cell membrane, ERa as positive 
control. Data are mean+S.E.M. from 3 independent experiments performed in 
duplicate (n= 3). 

(B-C) Calcium assay and IP1 assay to validate Lys05, but not estrogen related 
compounds, induces GPR30 activation and Gq signaling. Data are mean + S.E.M. from 


3 independent experiments (n = 3). 


Lys05-GPR30-Gq kDa 9,417 _Motion correction, 2D classification, 1,457,109 


movies Auto-picking (2 rounds) 
400 : 100 GPR30-Gq particles 
N 70 Ab-initio Reconstruction 2 round 
300 | Heterogeneous Refinement (2 rounds) 
~ alc: an 
é $ rih 
zA ORY ae 


274,405 particles 
Ab-initio Reconstruction 
(ese Refinement 
K uniform Refinement E 


2.5 
3.0 
3.5 
4.0 
5.0 


209,371 particles, 2.6 À 


ot 
TM7, e” 
TM1 GA 
Hah ) pe 
CAN ) 
= pe 


q 


UV280(mAU) 


5 10 15 


Volume(mL) 
(Superdex 200) 


H8 — Lys05-GPR30 — apo-GPR30 


Fig. S5 Purification and structure determination of Lys05-GPR30-Gq complexes. 
(A) Representative size exclusion chromatography (SEC) profiles and SDS-PAGE 
analysis of Lys-05 activated GPR30-Gq complex. 

(B-C) Representative cryo-EM image (B) and 2D classification averages (C) of Lys05- 
GPR30-Gq complexes. 

(D-E) Cryo-EM data processing flowcharts (D) and the Fourier shell correlation (FSC) 
curves (E) of Lys05-GPR30-Gq complexes. The global resolution of the final processed 
density map estimated at the FSC = 0.143 is 2.6 A. 

(F) The density maps of TM1—TM7 and helix 8(H8) of GPR30 in Lys05-GPR30-Gq 
complexes. 


(G) Comparison of Lys05 bound and apo-GPR30 structures. 
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Fig. S6 Structural analysis reveals a unique distribution of the transmembrane 
region of GPR30. 

(A) Structure comparison of GPR30’s 7TM bundle with that of SHT14 (PDB code: 
7E2Y), GPBAR (PDB code: 7CFM) and GPR183(PDB code: 7TUZ). 

(A) Structure comparison of GPR30’s 7TM bundle with that of GHSR (PDB code: 
7F9Z), NMUR2 (PDB code: 7W55) and OXR2 (PDB code: 7L1U). 
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Fig. S7 Mutational effects of Lys05-GPR30 interacting residues to the activation 
of GPR30. 

(A) Structure comparison of apo-GPR30 and Lys05-GPR30 indicated the 
conformational changes of related residues during Lys05 insertion. 

(B-E) Dose-response curves of Lys05 in activating the mutated GPR30, mutations are 
located in TM2 (B), TM3 (C), TM6 (D), and TM7 (E). Data are mean + S.E.M. from 3 


independent experiments (n = 3). 


Table S1 Cryo-EM data collection, refinement and validation statistics 


Voltage (kV) 

Electron exposure (e-/A2) 
Defocus range (um) 

Pixel size (A) 

Symmetry imposed 

Filtered particle images (no.) 
Final particle images (no.) 
Map resolution (A) 

FSC threshold 


Map resolution range (A) 


Refinement 
Initial model used (PDB code) 
Map sharpening B factor (A?) 
Model composition 
Non-hydrogen atoms 
Protein residues 
Ligands 
B factors (A?) 
Protein 


Ligand 


R.m.s. deviations 
Bond lengths (A) 
Bond angles (°) 

Validation 
MolProbity score 
Clashscore 
Poor rotamers (%) 
Ramachandran plot 
Favored (%) 
Allowed (%) 


Disallowed (%) 


GPR30-Gq complex 


E2 


300 
50 
-1.0 to -2.0 
1.071 
C1 
1,199,129 
177,591 
3.2 
0.143 
2.5 - 5.0 


AlphaFold GPR30 


G1 


300 
50 
-1.0 to -2.0 
1.071 
C1 
1,065,934 
411,141 
3.1 
0.143 
2.5-5.0 


AlphaFold GPR30 
-147.6 


Ful 


300 
50 
-1.0 to -2.0 
1.071 
C1 
3,410,733 
536,691 
3.2 
0.143 
2.5 -5.0 


AlphaFold GPR30 
-135.5 


apo 


300 
50 
-1.0 to -2.0 
0.824 
C1 
1,072,804 
300,561 
2.9 
0.143 
2.5-5.0 


AlphaFold GPR30 
-104.4 


Lys05 


300 
50 
-1.0 to -2.0 
0.732 
C1 
1,457,109 
209,371 
26 
0.143 
2.5-5.0 


AlphaFold GPR30 
-96.2 


Table S2. Effects of Lys05 on activation of GPR30 and its mutants. 


Receptor pECso ° Emax (%)* Surface Expression (WT%) ° 
WT 5.88+0.12 100+7.3 100+17 
L108 4.33+0.19 46.6+4.75 105+1.5 
Elis A UD* 10.5+7.0 142425 
L137 74 UD 13.3+1.0 143413 
Q138?%33 A UD 18.7+1.7 166+15 
M1412%°6A UD 12.7+6.8 114+5.1 
W272 A UD 8.03+41.5 89.4+4.8 
E2756%51A 4.9+0.07 80.6+8.0 1042.2 
H282°%8A 5.95+0.44 1042.2 1091.8 
H3007? A 5.68+0.21 40.0+7.1 81.7411 
P3037%7A 5.82+0.06 75.74+8.5 90.2+1.6 
N3107%*9A 5.75+0.14 84.445.5 98.2412 
F3147%8A UD 0.61+3.9 98.4+2.9 


è Data shown are means + S.E.M. from at least three independent experiments. 
> Data were normalized to the expression level of WT GPR30 in HEK293 cells. Data 
shown are means + S.E.M. from four independent experiments. 


“UD-undetectable, indicates that the activation level is too low to determine ECs values. 


